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Variability  of  the  Lower  Thermosphere  Determined 
From  Satellite  Accelerometer  Data 


i.  INTRODUCTION 

Accurate  density  measurements  have  been  previously  made  with  electrostatic 

12* 

accelerometers  on  low  altitude  satellites.  * These  in  situ  data  permitted  detailed 
analysis  of  atmospheric  behavior  compared  to  the  orbital  decay  technique  which 
involves  considerable  temporal  and  spatial  averaging.  Because  these  satellite  life- 
times were  relatively  short,  however,  only  limited  amounts  of  the  accelerometer 
data  were  available.  This  report  describes  an  extensive  neutral  density  data  base 
generated  using  data  from  accelerometers  on  four  low  altitude  satellites  (Atmosphere 
Explorer  -C,  -D,  and  -E,  hereafter  AE-C,  -D,  and  E,  and  Air  Force  satellite 
S3-1).  Data  obtained  from  over  4000  orbits  (a  total  of  115,951  discrete  values)  at 
altitudes  from  230  km  down  to  as  low  as  140  km  during  the  period  January  1974  to 
November  1976,  have  been  utilized.  This  represents  the  most  extensive  set  of 
neutral  density  measurements  in  existence. 


(Received  for  publication  25  May  1978) 

1.  Devries,  L.  L.  (1971)  Experimental  Evidence  in  Support  of  Joule  Heating 
Associated  With  Geomagnetic  Activity,  Marshall  Space  Flight  Center, 

S’ AS  A TVrX-GTSSBr 


Knowledge  of  'lospheric  density  and  its  variations  is  required  for  satisfactory 
low  altitude  satellite  design,  operation,  and  orbit  prediction.  Existing  models  of 
the  lower  thermosphere  are  deficient  since  they  are  primarily  based  on  data  obtained 
at  higher  altitudes.  This  report  compares  the  accelerometer  density  values  to  two 
commonly  used  models,  that  of  Jacchia1  hereafter  J71,  and  that  of  Hedin  et  al,  4,5 
hereafter  MS1S.  Comparisons  are  made  as  a function  of  altitude,  latitude,  geomag- 
netic activity,  and  local  time.  Frequency  distributions  of  these  comparisons  are 
described  in  terms  of  their  mean  value,  standard  deviation,  coefficient  of  skewness 
and  kurtosis.  Results  of  this  study  provide  a quantitative  estimate  of  the  improve- 
ments required  in  current  models  of  the  thermosphere.  Correlations  of  the  observed 
variations  with  solar  and  geophysical  parameters  will  be  described  in  forthcoming 
reports. 


2.  EXPERIMENT  AND  DATA  DESCRIPTION 


2.1  Satellite  Characteristics 

The  Atmosphere  Explorer  program  involved  three  low-altitude  satellites  de- 
signed to  permit  a coordinated  study  of  the  thermosphere.  Dalgarno  et  al^  have 
described  the  AE  mission.  AE-C,  -D,  and  -E  were  launched  into  low-perigee  high- 
eccentricity  orbits  with  different  inclinations.  A propulsion  system  was  used  to 
maintain  an  elliptical  orbit  for  approximately  1 year  (AE-D  data  acquisitions  were 
terminated  after  3-1/2  months  when  the  spacecraft  power  system  failed.  However, 
data  were  obtained  on  81  percent  of  the  orbits).  AE-C  and  AE-E  were  put  into  cir- 
cular orbits  near  250  km  following  their  elliptical  orbit  phase.  Data  in  this  report 
are  for  the  elliptical  orbit  phase  only.  Data  were  acquired  in  both  a spinning  (4  rpm) 
and  a despun  mode.  S3-1  was  also  launched  into  a low-perigee  high-eccentricity 
orbit  and  was  spin-stabilized  at  5 rpm.  This  satellite  had  a lifetime  of  nearly  seven 
months.  Table  1 provides  orbital  characteristics  and  the  dates  of  elliptical  orbit 


3.  Jacchia,  L.  G.  (197  1)  Revised  Static  -Models  of  the  Thermosphere  and  Exosphere 

With  Empirical  Temperature  Profiles.  Spec.  Rept.  332,  Smithsonian 
Astrophys.  Observatory,  Cambridge,  MA. 

4.  Hedin,  A.  E.  . Salah,  J.  E.  , Evans,  J.  V.  , Reber,  C.  A.  , Newton,  G.  P.  , 

Spencer,  N.  W.  , Kayser,  D.  C.  , Alcayde,  D.  , Bayer,  P.  , Cogger,  L.  , and 
.McClure,  J.  P.  (1977a)  A global  thermospheric  model  based  on  mass 
spectrometer  and  incoherent  scatter  data,  MSIS  1,  N„  density  and  temperature. 
J.  Geophys.  Res.  82t2139.  i 

5.  Hedin,  A.E. , Reber,  C.A.,  Newton,  G.  P.  , Spencer,  N.  W.  , Brinton,  H.  C.  , 

and  Mayr,  H.  G.  (1977b)  A global  thermospheric  model  based  on  mass 
spectrometer  and  incoherent  scatter  data,  .MSIS  2,  Composition,  J.  Geophys. 

Res.  82:2148.  

Vvw 

6.  Dalgarno,  A.,  Hanson,  W.  B.  . Spencer,  N.  W.  , and  Schmerling,  E.  R.  (1973) 

The  Atmosphere  Explorer  mission.  Radio  Sci.  8:263. 
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data  acquisition  for  each  satellite.  The  ait  tudes  at  which  data  were  obtained  ire 
shown  (shaded  areas)  as  a function  of  time  in  Figure  1.  The  AE-C  and  -E  perigees 
were  occasionally  lowered  to  below  140  km.  Perigee  of  S3-1  was  generally  between 
160  and  180  km.  Data  obtained  from  250  km  down  to  140  km  are  included  in  the 
data  base.  Figures  2a  through  2h  show  the  variation  in  latitude  (top  curve)  and  local 
time  (bottom  curve)  for  each  satellite. 

Table  1.  Satellite  Orbital  Characteristics  and  Data  Acquisition  Periods 


Satellite 

Launch 

Data 

End  Elliptical 
Data  Acquisition 

Inclination 

Initial 

Perigee 

Initial 

Apogee 

AE-C 

Dec  73 

Nov  74 

68° 

156  km 

4000  km 

S3- 1 

Oct  74 

May  75 

97° 

160  km 

4000  km 

AE-D 

Oct  7 5 

Jan  76 

to 

o 

o 

156  km 

3800  km 

AE-E 

Nov  75 

Nov  76 

o 

o 

CM 

157  km 

3000  km 
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1000  2000  3000  4000 


Figure  2a.  Latitude  and  Local 
Time  Distribution  for  AE-C 
160  km  Data 
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Figure  2b.  Latitude  and  Local 
Time  Distribution  for  AE-C 
240  km  Data 
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Figure  2c.  Latitude  and  Local 
Time  Distribution  for  AE-D 
160  km  Data 


RE-0  240  KM 


ORBIT  NO. 

400  800  1200 


ORBIT  NO. 

400  600  1200 


Figure  2d.  Latitude  and  Local 
Time  Distribution  for  AE-D 
240  km  Data 
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Figure  2e.  Latitude  and  Local 
Time  Distribution  for  AFJ-E 
160  km  Data 
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F igure  2g.  Latitude  and  Local 
Time  Distribution  for  S3-1 
160  km  Data 
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2.2  Drag  Measurement 

The  accelerometer  experiment  provides  extremely  accurate  measurements  of 
orbital  accelerations.  Three  single-axis  instruments  mounted  orthogonally  were 
flown  on  each  AE  satellite.  Operating  principles  of  this  experiment  have  been  given 

7 

by  Champion  and  Marcos.  S3-1  carried  a single-axis  sensor  essentially  identical 
to  the  AE  instrument.  The  sensitive  axis  was  aligned  to  the  flight  direction.  A 
one-quarter  sec  sample-time  was  used  for  each  instrument.  Drag  accelerations 

g 

were  separated  from  noise  accelerations  by  numerical  filtering  techniques.  For 
spinning  orbit  data  the  errors  are  estimated  as  follows:  area-to-mass  uncertainty 
is  ± 1%,  attitude  error  is  negligible  (only  ram  point  values  are  used),  and  filtering 
error  varies  from  negligible  at  perigee  to  ± 2%  at  200  km  to  ± 5%  at  250  km.  For 
despun  orbit  data  (AE  only)  the  sensors  are  at  45°  ± 2°  to  the  velocity  vector  and 
the  attitude  uncertainty  is  ± 2 %.  For  both  spinning  and  despun  data  there  may  be 
a ± 10%  error  in  the  assumed  free  molecular  flow  drag  coefficient  value  of  2.  2. 

The  drag  coefficient  error  and  area-to-mass  uncertainty  constitute  systematic 
errors.  Hence  the  random  error  in  measurement  of  density  variations  on  a spinning 
satellite  varies  from  negligible  at  perigee  to  ± 5%  at  250  km.  For  the  AE  despun 
mode  data,  due  to  the  accelerometer  alignment  with  respect  to  the  velocity  vector 
and  the  accuracy  of  the  attitude  control  system  there  is  an  additional  ± 2%  error. 
This  could  be  alleviated  by  aligning  the  accelerometer  sensitive  axis  with  the 
velocity  vector  or  by  a higher  accuracy  attitude  determination  system. 

7.  Champion,  K.  S.  W. , and  Marcos,  F.  A.  (1973)  The  triaxial  accelerometer 

system  on  Atmosphere  Explorer,  Radio  Sci.  8^:263. 

8.  Noonan,  J.P.,  Fioretti,  R.  W. , and  Hass,  B.  (1975)  Digital  Filtering  Analysis 

Applied  to  the  Atmosphere  Explorer-C  Satellite  ME&A  Accelerometer  Data, 


w oral  ta  M* 


2.3  Data  Ha*' 


Density  values  calculated  at  5-km  intervals  were  stored  in  the  data  base. 

Associated  with  each  point  is  the  corresponding  J7  1 and  A)  SIS  model  value. 

ephemeris  data,  solar  flux,  and  K value.  A histogram  of  the  data  distribution  as 

P 

a function  of  geographic  latitude  for  each  satellite  is  given  in  Figures  3a  through  3d. 
These  figures  show  the  number  of  data  points  obtained  in  each  10°  latitude  band. 

The  percentages  of  the  total  amount  of  data  (combining  the  four  sets)  in  the  latitude 
bands  0±  20°,  |20-50°|  and  |70-90°|  are  47%,  25%,  and  11%  respectively.  The 
large  amount  of  near  equatorial  data  is  due  to  the  low  inclination  of  AF-F  and  the 
precession  of  perigee  for  the  other  satellites  as  indicated  by  Figure  1.  The  data 
above  | 7 0 ° | are  provided  by  AF-D  and  S3-1.  Although  relatively  small  in  per- 
centage of  total  data,  12,822  points  were  obtained  at  these  latitudes.  Figures  4a 
through  4d  show  the  data  distribution  as  a function  of  local  time,  giving  the  number 
of  points  within  each  1-hr  interval.  The  figures  show  excellent  local  time 
coverage  provided  by  AF-C  and  -E.  The  distribution  of  data  with  respect  to  geo- 
magnetic activity  is  shown  in  Figures  5a  through  5d.  These  plots  give  the  number 
of  measurements  during  which  the  3-hourly  index  (with  a 6-hr  lag)  %11 
within  specific  1-unit  intervals.  Approximately  20%  of  all  the  data  occur  during 

conditions  when  the  average  K is  greater  or  equal  to  4 . The  relative  frequency 

P u 

distribution  for  K values  is  similar  to  that  obtained  from  a study  for  the  period 
9 P 

1932-1971.  In  Figures  6a  through  6d  distribution  of  data  with  solar  flux  is  shown. 
These  plots  give  the  number  of  measurements  during  which  the  Fjq  ^ value  (with 
a 1-day  lag)  fell  within  specific  five  solar  flux  unit  intervals.  It  can  be  seen  that 
the  data  base  was  obtained  mainly  during  conditions  of  very  low  solar  flux.  Only 


about  4%  of  the  data  occur  on  days  with  F 


> 100.  An  example!  of  the  data  base 


is  given  by  showing  results  obtained  at  the  altitude  of  180  km  for  each  satellite. 
F igures  7a  through  7d  show  these  density  data,  the  latitude  at  which  each  data 
point  was  obtained,  the  ratio  of  measured  density  to  the  J71  model,  and  the  K 


9.  C age,  A.  L.  , and  Zawalick,  F.  J.  (197  2)  A Discussion  of  the  Geomagnetic  Indices 
Kn  and  An,  1932  to  1971,  AFCRL-72-(J533: 
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RESULTS 


3.1  Statistical  Properties 

The  measured  density  values  have  been  statistically  analyzed  in  relation  to  the 
.171  and  MSIS  models.  The  normal  distribution  is  completely  determined  by  two 
constants:  the  mean  locates  the  center  of  the  distribution,  and  the  standard  devia- 
tion measures  the  spread  or  variation  of  the  individual  measurements.  Extreme 
deviations  from  models  are  of  particular  interest  to  potential  Air  Force  users, 
hence  departures  of  frequency  distributions  of  the  data  from  a normal  distribution 
have  also  been  calculated.  The  following  moments  of  x,  which  represent  the  per- 
cent departure  of  the  data  from  the  model  values  have  been  determined:^ 


Mean  value  (x) 


x * 


M 


£ 
i=  1 


(1) 


where  x^  is  the  percent  deviation  of  the  i*^1  density  measurement  from  the  model  and 
N is  the  total  number  of  data  points. 


Standard  deviation  (S) 


S = 


N 

£ 
i*  1 


(x.-x)2 

N-l 


1/2 


(2) 


For  a normal  distribution  about  68.3%  of  the  data  fall  in  the  interval  x ± S,  95.  5% 
are  within  x ± 2S,  and  99.73%  are  within  x ± 3S, 


Skewness  or. 

N 


i 


£ 

i=  1 


third  moment  about  the  mean, 
(x.-x)3 

*7  • 


(3) 


The  skewness,  a measure  of  non-symmetry,  will  be  positive  if  the  lower  values 

of  x are  closer  to  the  mean  than  the  higher  values,  that  is,  if  the  distribution  has 

a large  tail  to  the  right  (positive  side)  of  the  mean.  In  this  case,  large  positive 
3 — 

contributions  (x-x)  , when  x exceeds  x,  will  dominate  the  smaller  negative  con- 
tributions obtained  when  x is  less  than  x.  Negative  skewness  indicates  that  the 

lower  tail  is  the  extended  one.  The  factor  -4-  renders  this  measure  independent 

SJ 

10.  Snedecor,  G.  W.  , and  Cochran,  W.  G.  (1967)  Statistical  Methods,  Iowa  State 
University  Press,  Iowa. 
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of  scale.  If  yb~  is  zero,  the  distribution  is  symmetric  with  respect  to  the  mean. 
The  criterion  for  95%  confidence  that  the  non-zero  skewness  observed  is  not  due 
to  sampling  error  is  |^"b”.  | >(1.  where  N is  the  number  of  sampling  points 


Kurtosis  or  fourth  moment  about  the  mean,  (b2> 


7 


N 

s 

i=  1 


(x.-x)4 


(4) 


For  a normal  distribution  b,,  has  a value  of  3.  Values  of  b^  greater  than  3 
indicate  an  excess  of  values  near  the  mean  and  far  from  the  mean  with  a corres- 
ponding depletion  of  intermediate  values  of  the  distribution  curve.  For  values  of 
b2  < 3 the  distribution  curve  has  a flatter  top  than  normal.  The  criterion  for  95% 
confidence  that  the  kurtosis  observed  is  in  fact  different  from  3 and  not  due  to 
sampling  error  is  that 

|b2  - 3|><3.92)^. 

3.2  Frequency  and  Probability  Distribution 

Graphical  representations  of  frequency  distributions  of  percent  differences 
from  models  have  been  developed.  These  histograms  are  shown  for  each  set  of 
satellite  data.  This  separation  by  satellite  permits  some  intercomparison  as  a 
function  of  the  different  geophysical  conditions  appropriate  to  each  set  of  data. 

Planned  studies  of  atmospheric  variability  related  to  geomagnetic  activity,  latitude, 
longitude,  season,  and  the  semi-annual  variation  will  use  the  data  base  in  a unified 
format. 

Histograms  are  shown  in  F igures  8a  through  8d  for  each  set  of  satellite  data 
(at  all  altitudes)  compared  to  the  J71  model.  The  plots  show  the  percentage  of  the 
total  number  of  data  points  falling  within  each  1%  model  deviation  interval. 

The  total  number  of  points  is  given  at  the  bottom  of  the  box  in  the  upper  right  hand 
corner  of  each  figure.  A normal  distribution  curve  is  superimposed  as  a solid  curve 
on  each  histogram.  Deviations  from  the  normal  curve  are  observed  with  each  set 
of  satellite  data.  Values  for  the  coefficients  in  Eqs.  (1)  to  (4)  are  provided  in  the 
right  hand  corner  of  each  figure.  The  mean  values  fall  within  the  range  6.  9 t 2.  5%. 
The  standard  deviations  are  within  the  range  16.  1 ± 0.6%.  The  skewness  coefficient 
is  positive  for  all  cases  indicating  an  excess  of  large  positive  deviations  relative  to 
a normal  distribution.  The  kurtosis  has  a value  greater  than  3 for  each  case.  These 
properties  of  the  third  and  fourth  moments  are  statistically  confirmed  since 
1 . 96y^” and  3.92  a re  less  than  JT,  and  | bg- 3 1 respectively  for  each  case. 


25 


Figures  9a  through  9d  show  the  same  data  as  in  Figures  8a  through  8d,  but 
compared  to  the  MS1S  model.  This  model  gives  better  agreement  to  the  average- 
measured  absolute  density  values  than  .171.  The  mean  falls  within  2.9  ± 9.  6%  for 
each  set  of  data.  Standard  deviations  are  also  slightly  lower,  falling  within  the 
range  14.  1 ± 1.7%.  As  with  the  data  in  Figures  8a  through  8d,  the  newness  coef- 
ficient is  always  positive  and  the  kurtosis  is  always  greater  than  2.  No  trend  was 
found  for  these  coefficients  in  comparing  the  data  in  f igures  8a  through  8d  and  9a 
through  9d. 

The  MSIS  model  predicts  density  values  which  are  on  the  average  4%  lower  than 
those  of  the  .171  models.  Since  the  drag  coefficient  may  differ  systematically  from 
its  assumed  value  of  2.  2 by  about  10%,  it  cannot  be  determined  which  model  more 
accurately  predicts  the  absolute  density. 

An  earlier  study  comparing  OV1-15  satellite  accelerometer  data  to  .17  1 found 
an  average  departure  of  about  -1%  from  the  model.  These  data  were  obtained  for 
an  average  solar  flux  of  145  units.  However,  it  cannot  be  determined  whether  the 
.171  model  has  a systematic  error  with  respect  to  solar  flux  because  of  drag  coef- 
ficient uncertainties.  OV1-15  had  an  1,/I)^2  compared  to  AE  which  had  an  1,/D-l. 
For  the  higher  L/U  the  theoretical  value  of  ( ^ is  higher  by  about  10%. 

Cumulative  frequency  distribution  curves  are  another  useful  means  of  displaying 
the  deviations  of  the  measured  data  from  a normal  distribution.  Figures  10a  through 
1 0d  show  the  results  of  f igures  8a  through  8d  plotted  in  this  alternate  format.  The 
ordinate  of  a point  on  the  curve  gives  the  percent  of  the  total  data  which  deviates 
from  the  model  by  a percent  less  than  or  equal  to  the  value  of  the  abscissa.  The 
ordinate  scale  at  the  right  gives  the  cumulative  percent  in  terms  of  standard 
deviations.  These  curves  can  be  used  to  estimate  the  amount  of  expected  large 
* deviations  from  a model.  For  example,  approximately  2.  3%  of  all  the  AE-C'  data 

differ  from  the  J 7 1 model  by  between  -80%  and  -20%  and  approximately  2.  3%  of  this 
same  data  differ  from  the  model  by  more  than  45%.  The  normal  distribution  shows 
a linear  dependence  in  this  type  plot  and  is  shown  as  a solid  line. 

- 

i 

11.  Marcos,  F.A.,  and  Champion,  K.S.  W.  (1912)  Variations  of  the  neutral  atmos- 

pheric density  at  low  satellite  altitudes,  Free,  International  Conf.  fin  Aero- 
space and  Aeronautical  Meteorology , W ashington,  t).  C. 

12.  Karr.  G.A.,  and  Smith,  K.  H.  (1972)  Influence  of  satellite  aerodynamics  on 

atmospheric  density  determination,  Froc.  International  Conf.  on  Aerospace 
and  Aeronautical  Meteorology,  Washington,  IX  C. 
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3.3  l)ala  as  a Function  of  (leonuf'net it  Activity  ami  Latitude 

Frequency  distributions  were  made  for  data  obtained  at  different  levels  of  geo- 
magnetic activity  and  at  different  geographic  latitudes.  The  S3  - 1 data  were  used 

since  it  has  the  best  distribution  with  respect  to  K and  latitude.  The  data  are  first 

P 

displayed  as  a function  of  geomagnetic  activity  only.  Three  groups  of  data  were 
formed  according  to  the  index  (with  a 6-hr  lag).  Group  A included  measurements 
obtained  with  the  condition  s 3..  Groups  B and  C were  for  data  with  5 4(j 
and  Kji  5(  respectively.  The  range  4.  3*  K > 3-  was  omitted  to  provide  a larger 
difference  between  Groups  B and  C.  .Next  the  data  were  separated  into  30°  latitude 
bins.  The  mean  and  moments  about  the  mean  were  calculated  for  Group  A and  B. 

Results  obtained  when  the  data  are  separated  by  K only  are  shown  in  Figures 
1 la  through  11c  and  12a  through  12c.  These  figures  show’  the  frequency  distributions 
for  data  compared  to  J7  1 and  MSIS  respectively.  The  results  are  summarized  as 
follows: 

(1)  Mean  value:  The  actual  average  measured  increase  in  density  for  Group  B 
relative  to  Group  A is  18%.  l or  Group  C relative  to  Group  A it  is  24%.  The  J7  1 
Group  B and  C ratios  are  about  1.5%  less  than  those  of  Group  A.  The  MSIS  Group 
B and  C ratios  are  about  5%  greater  than  those  of  Group  A.  Hence  the  .171  model 
provides  a more  accurate  representation  of  the  geomagnetic  activity  effect.  The 
MSIS  model  underestimates  it  by  about  5%  probably  because  this  model  uses  a 
daily  average  value  to  estimate  the  density  enhancement  related  to  geomagnetic 
activity. 

(2)  Standard  deviation:  For  all  six  cases  S is  within  15.  3 ± 1%.  The  largest 
values  of  S occur  with  the  Group  C data  for  both  models. 

(3)  Skewness:  An  increase  in  bj  with  increasing  is  observed  in  Figures 
11a  through  11c  and  12a  through  12c.  This  is  indicative  of  a relative  increase  in 
the  number  of  positive  deviations  from  the  mean  with  increasing  K^. 

(4)  Kurtosis:  The  value  of  b„  decreases  as  K increases.  For  both  Group  A 

2 p 

sets  b^  = 5.  3.  For  Groups  B and  C 3.  4 s s 3.  9.  With  increasing  K the  peak  of 
the  distribution  flattens  out  while  the  number  of  intermediate  values  increases. 

The  Group  A and  Group  B data  were  next  separated  into  30°  latitude  bins.  Group 
C data  were  not  included  in  this  study  because  a relatively  small  number  of  points 
were  obtained  at  high  southern  latitudes.  Calculated  values  of  the  mean  and  mom- 
ents about  the  mean  have  been  tabulated  for  this  study  and  are  shown  in  Table  2. 

Table  2 shows  data  for  Groups  A and  B (low  and  high  K conditions,  respectively) 
relative  to  both  the  .17  1 and  MSIS  models,  for  low  conditions  (Group  A)  the  MSIS 
model  appears  to  provide  a more  accurate  depiction  of  the  latitudinal  structure;  x 
is  within  -4.  3 ± 2.  6%  for  MSIS  and  8.  3 ± 5.  9%  for  J7  1.  In  general,  J7  1 overesti- 
mates the  average  geomagnetic  activity  effect  by  about  4%  from  -60°  to  -30°  and  ac- 
curately depicts  it  from  +30°  to  -90°.  MSIS  values  are  about  6%  low  from -30°  to -60°  and 


32 


about  3%  low  in  the  latitude  intervals  -60°  to  +30°  and  + (>0°  to  +h0°.  The  largest 
standard  deviation  occurs  at  -30“  to  -00“  for  both  models.  However,  relatively  low 
values  occur  at  +30°  to  +60“.  As  was  found  in  f igures  11a  through  lie  and  I3a 
through  12c,  skewness  increases  with  increasing  K and  kurtosis  decreases. 


Table  2.  Mean  and  Standard  Deviation  as  Percent  Deviation  from  37  1 and  \1SIS 
Models  for  S3-1  Low  and  High  K ^ Data  Group  vs  Latitude 
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3.1  Altitude  Dependence 

Mean  and  standard  deviation  values  were  calculated  for  the  ratio  of  data-to- 
model,  and  the  average  percent  deviation  from  the  model  was  plotted  as  a function 
of  altitude.  Values  were  computed  at  5 km  intervals  using  all  data  for  each  satellite. 
Figures  13a  through  13d  and  14a  through  14d  show  results  compared  to  J7  1 and 
MSIS  respectively.  Mean  values  are  represented  by  crosses;  standard  deviations 
are  shown  as  solid  lines.  The  mean  values  with  respect  to  .17  1 are  between  - 6% 
(AE-D)  and  +15%  (S3 - 1 ) at  160  km.  Comparison  is  made  at  this  altitude  since 
below  160  km  the  number  of  data  points  decreases  and  varies  significantly  from 
satellite  to  satellite  (see  Figure  1).  At  240  km  x is  between  -2%  (S3  - 1 ) and  +4% 
(AE-E).  With  respect  to  MSIS  x is  between  ‘4%  (AE-U)  and  • 10%  (AE-C)  at  160  km 
and  -9%  (S3-1)  and  -2%  (AE-E)  at  240  km.  MSIS  provides  a particularly  accurate 
representation  of  the  AE-E  data.  The  .J7  1 scale  heights  are  in  general  10%  too  high 
below  about  170  km.  Above  this  altitude  they  are  about  10%  lower  than  the  measured 
data.  The  MSIS  scale  heights  are  also  about  10%  too  low.  Standard  deviations  are 
about  1 5%  at  160  km  and  increase  to  about  20%  at  240  km.  The  increase  is  due  to 
the  approximately  5%  filtering  error  in  the  data  reduction  at  this  altitude.  The 
large  standard  deviation  for  S3-1  occurs  below  160  km  where  there  are  relatively 
few  data  points. 


38 


* 


6E-C 


fif-0 


figure  l '.in.  Percent 
Deviation  from  the  .17  1 
Model  of  the  Mean  and 
Standard  Deviations  of 
AK-C  Density  .Measure- 
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Figure  13b.  Percent 
Deviation  from  the  .17  1 
Model  of  the  Mean  and 
Standard  Deviations  of 
AK-D  Density  Measure- 
ments vs  Altitude 
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figure  13c.  Percent 
Deviation  from  the  .17  1 
Model  of  the  Mean  and 
Standard  Deviations  of 
Af.-K  Density  Measure- 
ments vs  Altitude 
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Deviation  from  the  .17  1 
Model  of  the  Mean  and 
Standard  Deviations  of 
S3-  1 Density  Measure- 
ments vs  Altitude 
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figure  14a.  Percent 
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and  Standard  Deviations 
of  AK-C  Density  Measure- 
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figure  14h.  Percent 
Deviation  from  the 
MSIS  Model  of  the  Mean 
and  Standard  Deviations 
of  AE-D  Density  Measure- 
ments vs  Altitude 
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3. 5 Local  Time  Variation* 


A description  of  the  daily  variations  in  the  lower  thermospheric  density  was 
obtained  using  AE-E  data.  This  satellite  was  used  because  of  its  low  inclination 
and  excellent  local  time  coverage.  Data  at  160,  200,  and  240  km  were  separated 
into  1-hr  local  time  bins  and  averaged.  Similarly  an  average  ratio  to  .17  1 and  to 
MSIS  was  calculated  for  each  bin.  The  results  at  these  three  heights  are  plotted 
in  Figure  15.  Average  density  values  are  shown  as  solid  lines,  crosses  represent 
ratios  to  ,17  1,  and  circles  represent  ratios  to  MSIS.  Figure  15  (160  km)  shows  a 
density  minimum,  rather  than  maximum,  near  1400  hours.  Two  nearly  equal 
maxima  are  observed,  one  near  0700  and  one  near  2100  hours.  Evidence  of  the 
diurnal  bulge  is  seen,  increasing  with  altitude  at  200  and  240  km.  'These  results 
are  characteristic  of  atmospheric  tidal  variations.  ' * The  J7  1 model  does  not 
include  local  time  variations  due  to  tides.  The  density  is  underestimated  from  early 
evening  hours  to  late  morning  hours.  Ratios  vary  by  about  30%  for  the  data  shown 
in  Figure  15.  MSIS  does  include  tidal  components.  The  observed  variations  are 
more  closely  represented  and  the  agreement  is  within  +9  to  -1%,.  Figures  16a  and 
16b  show  an  example  of  the  frequency  distributions  for  a specific  local  time  interval. 
Data  obtained  at  all  altitudes  in  the  0800-0000  bin  are  shown.  Relative  to  .171,  ic  is 
11.2%  and  S is  11.9%.  For  MSIS  these  values  are  5.39%  and  11.92%,  respectively. 
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Figure  15.  Local  Time  Variation  of 
AE-E  Density  Data  at  160,  200,  and 
240  km;  Ratio  of  Data  to  .17  1 and 
MSIS 


13.  Mayr,  11.0.,  Harris,  I.,  and  Volland,  M. 
in  the  thermosphere,  ,1.  Oeophys.  Res. 


(1973)  Theory  of  the  phase  anomaly 
78:7480. 
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The  AK-K  data  are  further  examined  in  Figure  17  which  shows  the  density 
values  at  5 km  intervals  from  145-240  km  as  a function  of  local  time.  To  permit 

better  analysis  of  the  relative  tidal  amplitudes  the  data  have  been  normalized  to  an 

14 

altitude  of  ifiO  km.  Newton  et  al  have  detected  tidal  variations  in  the  altitude 
region  220  - 200  with  the  San  Marco  satellite  (inclination  3°).  Between  220  and  25 
km,  where  comparison  can  be  made,  their  results  agree  in  general  with  those  of 
Figure  16.  However,  the  semi-diurnal  component  appears  to  be  less  pronounced  in 
their  data.  The  difference  may  be  attributable  to  different  solar  and  geophysical 
conditions  for  AK-K  and  San  Marco.  The  AK-K  data  cover  a wider  latitudinal  and 
seasonal  range  and  were  obtained  during  lower  solar  flux  conditions  than  that  of 
Newton  et  al.  A report  describing  the  amplitude  and  phase  of  the  observed  AK-K 
local  time  variations  is  in  preparation. 

NORMALIZED 
i DENSItV 


Figure  17.  Altitude  Dependence  of  AK-K  I .oral  Time  Variations 
(Normalized  to  1(10  km) 


t.  koncm  sions 

An  extensive  data  base  has  been  developed  utilizing  accelerometer  measure- 
ments obtained  with  four  low  altitude  satellites.  In  addition  to  density  values  the 
data  base  incorporates  appropriate  satellite  orbital  information,  solar  and  geo- 
physical parameters  and  atmospheric  model  values.  A description  of  the  satellites 
the  accelerometer  experiment,  and  the  data  has  been  given.  It  is  anticipated  that 
this  data  base  will  he  useful  for  correlative  studies  of  aeronomical  problems  re- 
quiring knowledge  of  the  neutral  atmosphere.  The  present  study  has  been  directed 
toward  a statistical  analysis  of  atmospheric  variability.  Several  methods  of  dis- 
playing the  data  were  developed.  Analyses  included  determination  of  mean  values 

4.  Newton,  ( i.  1'.  , Kasprzak,  W.T.,  Curtis,  S.  A.  , and  I’elz,  D.T.  (16751  Local 
time  variation  of  equatorial  thermospheric  composition  determined  by  the 
San  Marco  3 Vice,  .1.  (ieophys.  Ites.  80:2280. 


and  the  second,  third,  and  fourth  moments  about  the  mean.  Comparison  has  been 
made  with  two  commonly  used  atmospheric  models,  .17  1 and  MSIS.  Results  have 
been  presented  in  Figures  8 through  17  and  Table  2.  These  data  give  the  extent  of 
unmodeled  variations  as  a function  of  geomagnetic  activity,  latitude,  altitude,  and 
local  time.  Neither  model  gives  a completely  adequate  description  of  the  observed 
atmospheric  variability  as  a function  of  all  these  parameters.  The  J7  1 model  gives 
a better  representation  of  the  average  geomagnetic  activity  effect.  MSIS  more 
accurately  describes  the  near  equatorial,  low  geomagnetic  activity  conditions.  Be- 
cause of  systematic  errors,  due  mainly  to  drag  coefficient  uncertainties,  it  cannot 
be  determined  whether  the  absolute  value  of  density  is  more  accurately  represented 
by  .171  or  MSIS.  Model  deficiencies  are  observed  in  the  depiction  of  local  time 

variations,  particularly  near  160  km  although  MSIS  provides  a better  depiction  than 

1 5 

.17  1.  An  evaluation  of  other  models  including  the  revised  .17  1 model  is  planned. 

The  data  base  will  be  utilized  to  develop  an  improved  model  of  the  lower  thermosphere 

1 6 

extending  the  results  of  Marcos  et  al  which  used  only  AE-C  data.  It  is  anticipated 
that  this  model  will  describe  the  density  data  in  this  report  as  a function  of  latitude, 
longitude,  altitude,  geomagnetic  activity,  and  local  time. 
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16.  Marcos,  F.  A.,  Garrett,  II.  B.  , Champion,  K.S.  W.  , and  Forbes,  .1.  M.  (1977) 
Density  variations  in  the  lower  thermosphere  from  analysis  of  the  AE-C 
accelerometer  measurements.  Planet.  Space  Sci.  25:499. 
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